Exact Nondipole Kramers-Henneberger Form of the Light-Atom Hamiltonian: An Application to 
Atomic Stabilization and Photoelectron Energy Spectra 
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The exact nondipole minimal-coupling Hamiltonian for an atom interacting with an explicitly time- and 
space-dependent laser field is transformed into the rest frame of a classical free electron in the laser field, i.e., 
into the Kramers-Henneberger frame. The new form of the Hamiltonian is used to study nondipole effects in 
the high-intensity, high-frequency regime. Fully three dimensional nondipole ab initio wavepacket calculations 
show that the ionization probability may decrease for increasing field strength. We identify a unique signature 
for the onset of this dynamical stabilization effect in the photoelectron spectrum. 



PACS numbers: 42.50.Hz,32.80.Rm,32.80.Fb. 

The general field of laser-matter interactions is character- 
ized by impressive progress in light-source technology. Light 
sources with pulses of shorter and shorter duration and ever 
increasing intensities are being developed. Pulses containing 
only a few cycles and with a duration of less than 1 fs are now 
commercially available Intensities of 10 14 W/cm 2 are 
routinely provided, and intensities two orders of magnitude 
higher, reaching the field strength of the Coulomb interaction 
in atoms and molecules, are not unusual. Femtosecond laser 
pulses have been used to produce coherent extreme-ultraviolet 
pulses of attosecond duration, and the expression "attosecond 
metrology" 1 2 j was coined for the investigation of matter with 
such short pulses |3]. Other developments include the large- 
scale intense free-electron laser projects at DESY (Hamburg, 
Germany) and SLAC (Stanford, USA). The TESLA test fa- 
cility in Hamburg has begun operation in the far-ultraviolet 
regime and, e.g., a study of the interaction of intense soft X- 
rays with atom clusters was reported ^j. The clusters ab- 
sorbed energy much more efficiently than anticipated from 
existing models, and the physical mechanism responsible for 
the excess in the absorbed energy is currently subject to some 
controversy @|. 

Typically the laser-atom interaction is described in the 
dipole approximation where several equivalent formulations 
exist; the most popular ones being the velocity gauge, the 
length gauge and the Kramers-Henneberger frame |6]. It is, 
however, clear that the new light sources alluded to above 
pave the way for studies of atomic and molecular systems 
under extreme nonperturbative conditions 0]. In the case of 
atoms interacting with light from the vacuum-ultra-violet free- 
electron laser the dipole approximation cannot be expected to 
be valid |8]. Thus, motivated by the need to include the full 
fc • r-term in the description of the light-matter interaction, we 
here revisit the question of equivalent formulations of electro- 
dynamics. 

We transform the exact nondipole minimal-coupling 
Hamiltonian for an atom in an explicitly time- and space- 
dependent field into the rest frame of a classical free elec- 
tron in the laser field. In the dipole approximation, this frame 
is known as the Kramers-Henneberger frame Our trans- 



formed exact nondipole Hamiltonian takes a simple form and 
is very useful for the discussion of strong-field dynamics. We 
apply it to the study of H in the high-intensity, high-frequency 
regime, and confirm the phenomenon of atomic stabilization, 
i.e., the possibility of having a decreasing ionization probabil- 
ity/rate with increasing intensity (for reviews see, e.g., I KID . 
Most importantly, we point out that the onset of the dynamic 
stabilization can be directly observed from electron energy 
spectra. [Atomic units (a.u.) with m e = e = h = 1 are used 
throughout. All derivations are straightforwardly generalized 
to atoms and molecules involving more electrons.] 

The minimal coupling scheme determines the Hamiltonian 
for a charged particle in an electromagnetic field through the 
vector potential A(rj) with rj = u)t — k ■ r, and k the wave 
number. The scheme implies that the canonical momentum is 
obtained by p — > p — qA and for an electron of charge q = 
— 1 in atomic units, we have p + A, and the time-dependent 
Schrodinger equation reads 

id t * v (r, t) = [(p + A{ v ) f + V(r)] *„(r, t), (1) 

where the subscript v refers to the velocity gauge. The ad- 
vantage of this formulation is that the spatial dependence of 
the field is explicitly accounted for through its presence in the 
vector potential. A disadvantage is that the interaction is not 
expressed in terms of the physical E- and S-fields. Also nu- 
merically, the evaluation of the action of the A ■ p-term can 
be quite involved unless a diagonal representation of $„ with 
respect to this operator is applied. Until now only the alterna- 
tive multipole formulation of Power-Zienau- Woolley 0, [Till 
has, in principle, kept the spatial dependence to all orders. 
The multipolar form represents the interaction in terms of the 
physical fields and the electron coordinate r, but, as the name 
suggests, it is inherently designed to provide an expansion of 
the light-matter interaction, and consequently very impracti- 
cal if one wishes to retain k ■ r to all orders. 

Here, we transform the Schrodinger equation into a new 
form by applying a nondipole Kramers-Henneberger transfor- 
mation. Let 

^kh = U^ V = cxp[i a(r)) ■ p\V v , (2) 
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FIG. 1 : Ionization and ground state probability for a two-dimensional 
model atom 1 16] in the nondipole (solid curve) and dipole (dashed 
curve) descriptions vs. electric field strength for a 5 cycle pulse with 
lo — 1 a.u.. 



p 2 + A 2 , s k 2 f da. 



which holds for a general elliptically polarized field. Within 
the dipole approximation A and a are space-independent, the 
last three terms are absent, and Q reduces to the well-known 
result 1 9 j . In the nondipole case, the importance of these terms 
is readily understood, e.g., in terms of their effect on a contin- 
uum wave function. The two terms proportional to k 2 are of 
the order of E 2 v 2 / (u> 2 c 2 ) and Eqv/c 2 , respectively, whereas 
the last term is of order Eqv 2 / (uc). We thus see that the ef- 
fect of the dominant term on a wave function is reduced by a 
factor ~ Eo/(ojc) compared to the p 2 -term. The factor Eq/lo 
is precisely the quiver velocity of the electron ^quiver, so we 
expect that the last three terms may be neglected as long as 
Vquiver/c <C 1. Whenever this condition is fulfilled, the non- 
relativistic approach is automatically justified as well. As it 
turns out, for the field parameters considered here, the effect 
of the nondipole terms is effectively given by the spatial de- 
pendence of the vector potential in the A 2 -term. 

As a first application of the new form of the Hamiltonian we 
consider the interaction with high-intensity, high-frequency 
fields. In this so-called stabilization regime [ 10], atoms may 
go through a region of decreasing ionization for increasing 
field strength. Stabilization was experimentally observed with 
Rydberg atoms 1 1 3 ] . With the development of new light 
sources, dynamic stabilization of ground state atoms is, how- 
ever, expected to be within experimental reach in the near fu- 
ture Q. 

Nondipole terms were investigated in approximate ways 
earlier and found to have a detrimental effect on the stabi- 



where 

a(r?) = - r drf A(rf) (3) 

represents the quiver motion relative to the laboratory frame 
of a classical free electron in the field. The Hamiltonian cor- 
responding to the new point of view is obtained by taking the 
time-derivative on both sides of @, and by using Q for ^> v , 
we obtain idt^ kh{j ', t) — Hkh^'kh(''', t) with 

H KH = UH v tf + i(d t U)tf. (4) 

To evaluate the effect of the unitary translation operators 
in (0}, we use the operator identity known as the Baker- 
Hausdorff lemma 111 211 and take advantage of the Coulomb 
gauge restriction [p, A] = and k ■ A = 0. The resulting 
Hamiltonian reads 



\ 2 ik 2 d 2 a ( da \ . 

■ p ) + - w - p+ UH (fe - p) ' (5) 
i 

lization lll5l[l6ll . The relative role of the different nondipole 
terms in Q is illustrated in Fig. ^ f° r a two-dimensional 
model atom 11611 . The ground state was exposed to a laser 
pulse propagating in the x direction and of linear polariza- 
tion u, p along the z-axis corresponding to the vector poten- 
tial A(r~i) = ^jf{rj) sin(?7 + <p)u p with f(ri) the envelope, 
Eq the electric field amplitude and <f> a phase that ensures 
that the vector potential represents a physical field 11711 . The 
wave function was propagated on a Cartesian grid by means 
of the split-step operator technique . A 5-cycle laser pulse 
with central frequency lo = 1 a.u. (46 nm) corresponding 
to the pulse duration T = 760 as, and with carrier-envelope 
f(rj) = sin 2 (^y), was employed. The intensity range was 
set to < I < 1.4 x 10 19 W/cm 2 . The population not ac- 
counted for in Fig.^ is left in excited states. The effects of the 
last three terms in the Hamiltonian as well as the spatial de- 
pendence of the quiver amplitude a(r, t) cannot be resolved 
on the scale of Fig.^ Hence, the nondipole effect observed 
for field strengths greater than 7 a.u. are exclusively related to 
the spatial dependence of the A 2 -term. 

We have, accordingly, justified that for the parameters un- 
der concern, it is a very accurate approach to apply the Hamil- 
tonian 10, neglecting the last three additional kinetic energy 
terms arising from the transformation (|2ji, to a fully three- 
dimensional study of ionization of a real ground state atom by 
intense short wave light field beyond the dipole approxima- 
tion. We consider H(ls) exposed to 5-cycle laser pulses in the 
attosecond range with central frequencies lo — 1 a.u. and ui = 
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FIG. 2: Upper panel: As Fig.Q but for the fully three-dimensional 
case with the system initially prepared in the H(ls) state. Lower 
panel: corresponding results for u = 2 a.u.. 

2 a.u.. The time-dependent Schrodinger equation is solved 
numerically based on a split-step operator approximation on a 
spherical grid as detailed elsewhere 1 19]. The wave function 
is expanded on the grid points [(r^iljfc) = (r^, 9j, <fik)] as 

ty(ri,ti jk ,t) = fi,m(n,t)Yi jm (Q jk ), (6) 

i. ni 

and the initial field-free H(ls) state is obtained from the exact 
analytical expression. Reflection at the edges r — r max — 
200 a.u. is avoided by imposing an absorbing boundary. For 
convergence, we include harmonics up to l ma x = 29, check 
for gauge invariance, use propagation time-step At = 0.01 
a.u., and set Ar = 0.2 a.u.. Photoelectron probability dis- 
tributions are calculated by projecting the wave function onto 
the field-free (discretized) continuum states. We note that the 
presence of nondipole terms will lead to a population of dif- 
ferent m-values in (jfj. 

In Fig.|2]total ionization- and ground state probabilities are 
shown for the fully three-dimensional case in the nondipole 
and dipole limits for two different frequencies. We ob- 
serve that the dipole approximation remains valid up to field 
strengths of the order of 10 a.u., and we find in general only a 
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FIG. 3: Upper panel: Ionization probability density dP/dE for 
H(ls) vs. electric field strength and energy of the ionized electron 
for a squared shaped 5-cycle pulse (380 as) with uj — 2 a.u. (23 
nm) by fully three-dimensional calculations. Middle panel: Results 
from the time-averaged potential 0. Lower panel: Total ionization 
probability in the dipole approximation (full curve) and for the time- 
averaged potential (dashed curve). 

small effect of the nondipole terms on stabilization. 

We now turn to the central question of how stabilization 
can most efficiently be experimentally detected. Measurement 
of absolute probabilities will require control of all parameters 
of the experiment: atom density, pulse characteristics, repe- 
tition rates, electron counts etc.. We therefore suggest mea- 
suring the energy-differential photoelectron spectrum. Figure 
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|3] shows the ionization probability density dP/dE vs. elec- 
tric field strength and energy E of the ionized electron with 
the full interaction potential (upper panel) and with the time- 
averaged Rramers-Henneberger potential (middle panel) IXOffl . 

V (a ;r) = i / V(r + a)dt, (7) 
1 Jo 

where a = E q /lj 2 is the quiver amplitude. The dipole and 
nondipole results are practically identical, and only the dipole 
results are shown in Fig. [5] For lower field strengths a reg- 
ular pattern of multiphoton resonances corresponding to ab- 
sorption of loo, 2ui or 3uj from the field is present. However, 
the multiphoton ionization process weakens at higher inten- 
sities as the stabilization sets in. Simultaneously, there is a 
steady growth in the portion of low-energy photoelectrons in 
the spectrum which can be assigned to Vq of 0. That Vq 
is responsible for the growth in the low-energy spectrum is 
readily seen by comparison of the upper and middle panels. 
The processes leading to ionization effectively divide into two 
competing classes: The multiphoton ionization superimposed 
on a monotonically increasing 'background' ionization pro- 
cess due to Vq solely. This is explicitly illustrated in the lower 
panel of Fig. [5] where the total ionization probability vs. elec- 
tric field strength is shown 1 20] . Multiphoton ionization dom- 
inates at lower field strengths, whereas the picture is the op- 
posite at higher values of Eq. The ionization due to the Vq 
potential reflects to what extent the laser pulse is turned- on 
and off adiabatically, and in a 'sudden approximation' picture 
it represents the lack of overlap between the field-free and the 
field-dressed states. Common in both photoelectron spectra 
is the presence of peaks in the probability density which can- 
not be attributed to multiples of w. Instead, they are a result 
of the non-adiabatic turn-on and turn-off of the field and can 
be associated with the higher-order Fourier components of the 
pulse. 

In summary, we presented a new formulation of the interac- 
tion between atoms and light maintaining full spatial depen- 
dence of the fields. We analyzed the terms in the interaction 
Hamiltonian and argued, supported by numerical evidence, 
that certain terms can be neglected. For the present field pa- 
rameters, the main nondipole effects come from A(rj) 2 . As an 
application, we considered the phenomenon of dynamic sta- 
bilization in intense high-frequency fields. We showed by full 
three-dimensional wavepacket simulations that the nondipole 
terms do not destroy the stabilization effect, and most impor- 
tantly that the photoelectron spectra in the stabilization regime 
shows very characteristic features: After onset of stabilization 
all ionized electrons have very low kinetic energy. Thus, by 
simply measuring the energy of the released electrons stabi- 
lization can be detected. 
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